LABORATORI NAZIONALI Dl FRASCATI 

SIS-Pubblicazioni 



LNF-04/07 



NEUTRINO OSCILLATION STUDIES WITH LASER-DRIVEN 

BEAM DUMP FACILITIES 

S.V. Bulanov 15 , T. Esirkepov 1,6 , P. Migliozzi 2 , F. Pegoraro 3 , T. Tajima 1 , F. Terranova 4 
1 - ) Kansai Research Establishment, JAERI, Japan. 
^INFN, Sez. di Napoli, Napoli, Italy 
^Dip. di Fisica, Univ. di Pisa and INFM, Pisa, Italy 
^INFN, Laboratori Nazionali di Frascati, Frascati, Italy 
5 ) A. M. Prokhorov Institute ofRAS, Moscow, Russia 
6 ) Moscow Institute of Physics and Technology, Dolgoprudny, Russia. 



Abstract 

A new mechanism is suggested for efficient proton acceleration in the GeV energy range; 
applications to non-conventional high intensity proton drivers and, hence, to low-energy 
(10-200 MeV) neutrino sources are discussed. In particular we investigate possible uses 
to explore subdominant — > v e oscillations at the atmospheric scale and their CP conju- 
gate. We emphasize the opportunity to develop these facilities in conjunction with projects 
for inertial confined nuclear fusion and neutron spallation sources. 



PACS: 14.60.Pq, 52.38.Kd, 41.75.Lx 



1 Introduction 



In the light of the strong evidence [1] for neutrino oscillations coming from atmospheric 
and solar neutrino experiments, recently corroborated by reactor and accelerator results, 
a very peculiar texture of the leptonic mixing matrix is emerging. Current results point 
towards two hierarchical mass scale differences 1 (Am^ <C (Am^l ~ jAm^D driving, 
respectively, the oscillations at the "solar" and "atmospheric" scales. The corresponding 
mixing angles are large but the interference between the two scales (sub-dominant — > 
v e oscillations at E/L ~ | Am| 3 |, L and E being the neutrino path-length and energy) has 
never been observed [2,3]. The leptonic mixing matrix (Pontecorvo-Maki-Nakagawa- 
Sakata, PMNS [4]) is usually parametrized [5] as: 



with Sij = sin 6ij and = cos 9^; current data suggest [6] at 90% confidence level 
35° < 9 23 < 55° , 0i2 = 32.5° ± 2.4° and small values for 13 (< 10°) i.e. support a 
"bi-large" PMNS. In fact, in the limit 13 — > the matrix becomes real- valued and the 
complex CP violating phase turns out to be unobservable. It follows that the possibil- 
ity to determine experimentally the (1, 3) sector of PMNS - i.e. the off-diagonal factor 
U e3 = s 13 e~ t5 - and, in particular, the Dirac phase 5 critically depends on the size of 9 13 . 
No theoretical inputs are available to constrain the size of 613 and S in a convincing man- 
ner, so that its experimental determination is mandatory. Such determination can be car- 
ried out at accelerators either measuring the size of the subdominant — > v e {p^ — > v e ) 
oscillation probability at the atmospheric scale or its T-conjugate v e — > (P e — > 9^). 
This measure will likely be the most challenging task of future long-baseline neutrino 
oscillation experiments and, for 613 values significantly smaller than current limits, tradi- 
tional accelerating techniques will be unable to provide the requested intensity and purity. 
Therefore, interest for unconventional neutrino sources has steadily grown in recent years, 
bringing e.g. to the proposals of the Neutrino Factories [7] or the Beta-beams [8]. 

Along this line, in this paper we consider a wide synergic scenario between nu- 
clear and neutrino physics programs connecting the long term development of facilities 
for laser-driven inertial confinement fusion (ICF) and the possibility to obtain an ultra- 
intense low-energy (1-2 GeV) proton driver for neutrino studies and spallation neutron 

1 We assume here only three active neutrinos with masses mi, to 2 and m 3 ; hence, only two independent 
mass scale differences exist (Amf 2 =m\—m\ and Am| 3 = m§ — m|) since Am^ = Am^ + Amf 2 - 



^(012,013,023,5) 
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sources (see Appendix). Such a link is made possible if an efficient laser-based acceler- 
ation mechanism is available in the near future. This mechanism is described in details 
in Sec. 2. The connection with neutrino physics results from the following considera- 
tion. Sources of neutrinos with energy beyond the /j, production threshold in charged- 
currents (CC) interactions are not strictly necessary to explore the magnitude and phase 
of U e3 . In principle, a high intensity source of and with energy of the order of a 
few tens of MeV would suffice to search for — > v e appearance and its CP-conjugate at 
baselines of L ~ 10 km. The neutrinos come from pion and muon decays at rest (DAR) 
and the pions can be produced by a high current proton beam dump facility. These beams, 
if available with proper intensity would provide simultaneously a source for — > v e os- 
cillations through the n + — > DAR chain and a source for P M — > v e through the 
subsequent ji + — > e + v^v e decay. Moreover, differently from Superbeams [9,10], the in- 
trinsic v e (P e ) beam contamination can be kept easily below 0.1% (see Sec. 3), as the 
facility can be operated below the K production threshold. 

The acceleration of protons up to several tens of MeV by the interaction of ultra- 
intense laser beams with solid targets has been recently reported by several experiments 
[11]. This process will probably open up a wealth of applications: radioisotope produc- 
tion [12], proton probing [13] and oncological hadron-therapy [14] have been discussed. 
Moreover, the laser-driven acceleration mechanism is a natural candidate for fast beam 
injection into conventional accelerators [15]. Part of these research programs might be 
carried out already with present laser technologies, provided that a suitable repetition 
rate, reproducibility and an improved beam quality become available. On the other hand, 
future developments towards higher laser intensities [16,17] would allow particle acceler- 
ation for High Energy Physics applications. 

In fact, laser-driven acceleration mechanisms mainly favor applications that profit 
of the significant beam intensity without imposing strong constraints on the beam quality 
(energy spread and emittance) and particle energy. The possibility of using the well- 
controlled time structure of laser-induced v beams to test the KARMEN time anomaly or 
increasing the laser intensity to overcome the kinematical threshold for muon production 
(y^ disappearance tests) has already been considered in [18]. The pion generation by the 
laser accelerated ions has also been discussed in [19]. Similarly, in this paper we suggest 
that the above-mentioned technique could be implemented to overcome the present lim- 
itations of proton dump facilities (Sec. 2) and, particularly, high intensity DAR neutrino 
beams (Sec. 3). We investigate its capability to clarify the (1,3) sector of the PMNS ma- 
trix, with emphasis on the size of 9 13 (Sec. 4-6), and the main technological challenges 
for a new generation of beam dump facilities. 
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2 Laser ion accelerators 

The classical mechanism of ion acceleration through the interaction of a laser pulse with 
a plasma target is a direct consequence of the electron acceleration. Due to the smaller 
electron mass, the energy of the laser light is first transformed into electron kinetic energy. 
The resulting displacement of the electrons and the modification of their density lead 
to the formation of a region with a strong electric charge separation. This causes an 
intense electric field which eventually accelerates the ions. In the simplest case of a 
one-dimensional geometry, when a transversally wide laser pulse interacts with a thin 
foil, the ponderomotive pressure of the laser pulse shifts the electrons with respect to 
ions, forming a strong electric field layer between the two species. An electric field also 
can be formed due to charge separation if the laser radiation accelerates a relatively large 
portion of the electrons, expelling them almost isotropically. In this case the fast electrons 
leave the targets and the heavier ions remain at rest forming extended regions of positive 
electric charge. Then the ions with non-compensated electric charge expand and acquire 
a kinetic energy. This corresponds to the so called "Coulomb explosion" regime. In 
configurations with more than one dimensions, different effects come into play such as the 
finite size of the waist of the laser beam, the transverse filamentation instability of a wide 
electromagnetic packet in a plasma, and the transverse modulations of the electron and 
ion layers. These effects usually reduce the energy of the fast ions and/or the efficiency 
of the energy transformation compared to the one obtained within the framework of the 
one-dimensional approximation. 

However, the process of ion acceleration exhibits new properties in the regime 
where the radiation pressure of the electromagnetic wave plays a dominant role in the 
laser-foil interaction, as demonstrated in Refs. [17,20]. In this regime electrons gain their 
energy due to the radiation pressure in a way that, qualitatively, resembles the ion accel- 
eration mechanism proposed by Veksler [21] in 1956. Here, the ion component moves 
forward with almost the same velocity as the average longitudinal velocity of the elec- 
tron component, hence with a kinetic energy well above that of the electron component. 
The ion acceleration appears to be due to the radiation pressure of the laser light on the 
electron component with the momentum being transferred to the ions through the electric 
field arising from charge separation. This mechanism of ion acceleration can be called the 
Radiation Pressure Dominant (RPD) mechanism. In contrast to the schemes previously 
discussed in the literature (see e.g. [22-24]) here the ion beam generation is highly effi- 
cient, and, as we will see later on, the ion energy per nucleon is proportional to the laser 
pulse energy. 
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Figure 1 : Interaction of the e.m. wave with the co-propagating proton-electron mirror. 

2.1 Radiation Pressure-dominated (RPD) regime: ID analytical description 

The acceleration mechanism can be explained as follows. The accelerated foil, which 
consists of the electron and proton layers, can be regarded as a relativistic plasma mirror 
co-propagating with the laser pulse. Assume that the laser pulse is perfectly reflected 
from this mirror. As a result of the reflection at the co-propagating relativistic mirror, the 
frequency of the electromagnetic wave decreases by the factor of (1 — v/c) / (1 + v/c) ~ 
1 / A'-/ 2 , where v is the mirror velocity and 7 = ( 1 — v 2 /c 2 ) ~ 1 ^ 2 is the Lorentz-factor of the 
plasma mirror 2 . 

Before the reflection, in the laboratory reference frame, the incident laser pulse 
energy £i as ,in is proportional to ElLi as ^ in , where E is the electric field amplitude and 
Lias,in is the incident pulse length. After the reflection the pulse energy becomes much 
lower: Bi as , re} oc Ef asref L las ^ reS « E$L lasM /4<y 2 . The length of the reflected pulse is 
longer by a factor A'j 2 , and the electric field is smaller by a factor 47 2 . Hence the plasma 
mirror acquires the energy (1 — l/47 2 )^ as in from the laser. In this stage the radiation 
pressure of the light accelerates the plasma slab (electrons and protons). As discussed 
above, the radiation momentum is transferred to the protons through the charge separation 
field and the kinetic energy of the protons is much greater than that of the electrons. This 
scenario is illustrated in Fig. 1, where a sketch of the e.m. wave interaction with the 
co-propagating proton-electron slab is presented. 

It is possible to estimate the proton maximum energy and the acceleration efficiency 

2 To ease comparison with existing literature, we use here the Gaussian CGS system of unit (c ^ 1). We 
resume the usual c = 1 system in the subsequent sections. 
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using the model of the flat foil driven by the e.m. radiation pressure, as described above. 
The radiation pressure is given by: 

2tt \u J 2tt 2tt \c + vj ' V 
where v = dx/dt is the foil instantaneous velocity, i.e. 

dx p 



dt (mlc 2 + p 2 ) 1 / 2 ' 



(3) 



Here primed (unprimed) quantities refer to the moving (laboratory) reference frame re- 
spectively. In a quasi-one-dimensional geometry, the laser electric field at the foil location 
x(t) depends on time as E' = E (t — x{t)/c). 

Since the expression of the radiation pressure is the same in both frames, we can 
write the equation of motion of the foil as 

dp = E 2 (t-x(t)/c) (m 2 p c 2 + p 2 ) 1 / 2 - p 

dt 2irn l (m 2 p c 2 + p 2 ) 1 / 2 +p' { 

where p is the momentum of the proton representing the foil, l and n are the thickness 
and initial proton density of the foil. In this approximation we neglect the heating of the 
proton fluid. 

In the simplest case, when the laser pulse is assumed to be long enough and with a 
homogeneous amplitude, we can consider the electric field E in Eq. (4) to be constant. 
Its solution p(t) is an algebraic function of time t. For the initial condition p = at t = 
it can be written in the implicit form 

2p 3 + 2(m 2 p c 2 + p 2 ) 3 ' 2 + 2 mc _ E l t (5) 
3m 2 c 2 ^ 3 p 27cn l 

At the initial stage, for t 2imolom p c/ E%, the proton momentum is a linear function of 
time: 

p w (E 2 /27rn / )t. (6) 

As t — > oo, the dependence of the accelerated proton momentum on time changes asymp- 
totically to 

p ~ £ p kin/c ~ m p c(3£ , Qt/87rno/o m p c ) 1/ ' 3 - (7) 

We notice here the obvious analogy between the proton motion regime described by the 
expression (5) and the solution of the problem of the acceleration of a charged particle 
under the radiation pressure of the electromagnetic wave. This analogy can be clearly 
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seen by comparing Eq. (5) with the expression for the velocity of an accelerated electron, 
which, in this limit, can be cast in the form Wa T t/m e c xs (2/3)(l — v 2 /c 2 ) 3 ^ 2 , with 
cW = cEq/Att the wave intensity (see Ref. [25]), a T = 87rrg/3 the Thomson cross 
section and r e = e 2 /m e c 2 the classical electron radius. In our case 2/n l plays the 
role of the effective cross section of the scattering of the electromagnetic wave. This 
analogy further underlines the similarity between the RPD mechanism of ion acceleration 
discussed here and the Veksler mechanism mentioned above. 

To find an upper limit to the proton energy acquired during the interaction with a 
laser pulse of finite duration, we must include the dependence of the laser electromagnetic 
field on time t and on the coordinate x. Because of the foil motion, the interaction time 
can be much longer that the laser pulse duration Tj aSjin . We introduce the phase of the 
wave 

ij> = u) (t-x(t)/c), (8) 

as a new variable, lu being the incoming laser frequency. Using Eq.(3) together with 
Eq.(8), we cast Eq.(4) for the particle momentum to the form 

dp _ E 2 {ip) (m 2 p c 2 +p 2 y/ 2 
dip 2-KLd n l (m 2 c 2 + p 2 ) 1 / 2 + p ' 

Its solution reads 

w(w + l) 

p = mpC ^TW) ' (10) 

where w is a function of ip: 

r a a*. (id 

Using Eqs.(8) and (10) we find the dependence of the foil coordinate on time and write 
the equations for t and x as functions of the variable ip in the form 



dt 1 f c \ 1 (m 2 p c 2 + p 2 ) 1 ' 2 



dip uj \c — v J uoq (m 2 c 2 + p 2 ) 1 / 2 — p 

= —[l + 2w(w + l)] (12) 

dx dx dt c „ , _ , j ^ 

dip dt dip LUq 

For a constant amplitude laser pulse with Eg(ip) = E^Oiip), where 9(ip) = for ip < 
and 9{ip) = 1 for ip > is the unit step function, we have w(ip) = w 9(ip)ip, with 



and 
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wq = El/ 1 (A-Kuoonolorripc). Then Eqs. (12) and (13) yield the parametric dependence of 
the accelerated foil coordinate on time 

t = (*/> + ttW> 2 + 2u> V73)M), (14) 

x = (w ij 2 + 2w 2 ij 3 /3)(c/u ). (15) 

In the limit t <C 2irn l m p c/ El we have x a cw t 2 /uj , and for t — > oo, x a ct, while 
the momentum p increases according to Eq. (7). 

The function w(tp) given by Eq. (11) can be interpreted as the normalized energy 
of the portion of the laser pulse that has interacted with the moving foil by time t. Its 
maximum value is w max = £i as ^ in /Af p m p c 2 , where 

E ScTl aSy i n 

C-las,in — : (tOJ 

47T 

is the laser pulse energy, Af p = n l S is the number of protons in the region, with area 
equal to S, of the foil irradiated by the laser pulse. 

From the solution of Eq. (4) given by Eq. (10) in terms of w we obtain for the 
kinetic energy of a proton initially at rest 

£ P kin = (rn 2 p c 2 + p 2 ) 1/2 c - m p c 2 = m p c 2 w 2 /(w + 1/2). (17) 

In the limits w < 1 and t« » 1, we have respectively 

£ P kin ~ 2m p c 2 w 2 , and £j,km ~ m p c 2 w. (18) 

The upper limit to the proton kinetic energy and, correspondingly, to the efficiency of 
the laser-to-proton energy transformation can be obtained from Eq. (17) by setting w = 

W max £las,in/ J^fpTflpC . 



2£las,in £las,in 



r> _ *"^las,in ^las,in /'1Q\ 

£-pkin,max 7w-> I T7 o T7 " 



We see that within this model almost all the energy of the laser pulse is transformed into 
the energy of the protons if £i as ,in 3> N p m p c 2 /2: 



-2 

J\f p 4:im lo 



F „ Slas,in „ E 2 CT las , ln 

^"pkin,max — K r — A , K.^ 1 ") 



It is worth noting that, when the RPD mechanism takes place, the dependence of the 
proton kinetic energy on the laser intensity and duration turns out to be linear. Moreover, 
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£ P kin,max does not depend on the size of the illuminated area S: within this simplified one- 
dimensional approximation, an increase of the focal spot S and the laser energy £i as , in 
that keep constant the ratio 

~ 4tt { } 

result in an increase of the number of accelerated protons J\f p = n l S without perturbing 
the proton energy spectrum. The acceleration length x acc ~ ct acc and the acceleration 
time t acc can be estimated using Eq.(20) and the t 1 / 3 asymptotic dependence of the proton 
energy on time in Eq.(14) as 

_ 2 / S lasM \ 2 

tacc ~ 7T I T7 o I 7~las,in- (22) 

3 \Al p m p c 2 J 

2.2 Computer Simulations 

Within the framework of the simplified ID approximation used above, the protons for- 
mally have a monoenergetic spectrum. A number of processes such as the transverse 
inhomogeneity of the amplitude of the laser pulse, the electron stochastization due to 
"vacuum heating" [26] and the subsequent proton layer expansion under the action of the 
Coulomb repelling force may result in the broadening of the proton energy spectrum or 
in inefficient acceleration. 

In order to examine this scheme in a three-dimensional geometry, whose effects can 
indeed play a crucial role in the dynamics and stability of the plasma layer under the action 
of a relativistically strong laser pulse, we have performed 3 3D particle-in-cell (PIC) sim- 
ulations with the code REMP (Relativistic Electro-Magnetic Particle-mesh code). This 
code is based on the current assignment scheme "Density decomposition" [27]. In these 
simulations the laser pulse is linearly polarized along the z-axis and propagates in the 
direction of the rr-axis. Its dimensionless amplitude is a = eE /m e ujc = 316, corre- 
sponding to the peak intensity I = 1.37 x 10 23 W/cm 2 x (1 /im/ A) 2 , A = 1/im being the 
laser wavelength. The laser pulse is almost Gaussian with FWHM size 8A x 25A x 25A 
and a sharp front starting from a = 100; its energy is £ L = 10 kJ x (A/1 /im) 2 . Protons are 
generated from a 1 /im foil. In fact, due to the finite contrast of the laser, pulse pedestals 
reaching the target before the main pulse will pre-form a fully ionized plasma [28,29]. 
Therefore, in the present situation the target behaves as a fully ionized, 1A thick plasma 

3 3D fully relativistic simulations of laser-plasma interactions represent major numerical efforts. The 
simulation under consideration has been performed on 720 processors of the supercomputer HP Alpha 
server SC ES40 at JAERI Kansai and form the basis of the studies of Refs. [17]. Here we make use of these 
results and, when appropriate, perform extrapolations to other flux and energy ranges. A more detailed 
numerical study is in progress and results will be presented in a forthcoming publication. 
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with density n e = 5.5 x 10 22 cm~ 3 x (1 /im/A) 2 , which corresponds to the Langmuir 
frequency uo pe = 7lu. The protons and the electrons have the same absolute charge and 
their mass ratio is m p /m e = 1836. The simulation box size is 100A x 72A x 72A corre- 
sponding to the grid size 2500 x 1800 x 1800, so the mesh size is 0.04A. The total number 
of quasi-particles is 4.37 x 10 9 . The boundary conditions are periodic along the y- and 
z-axis and absorbing along the x-axis for both the e.m. radiation and the quasi-particles. 
The results of these simulations are shown in Figs. 2-5, where the space and time units 
are respectively the wavelength A and period 2-k/lu of the incident radiation. 

Fig. 2 shows the proton density and the E z component of the electric field. We see 
that a region of the foil with the size of the laser focal spot is pushed forward. Although 
the plasma in the foil is overcritical, it is initially "transparent" for the laser pulse due to 
the effect of relativistic transparency (see e.g. [22]). Therefore a portion of the laser pulse 
passes through the foil. Eventually the pulse accelerates the electrons and, as a result of 
the charge separation, a longitudinal electric field is formed. This can be interpreted as 
the "rectification" of the laser light, by analogy with a rectifier in electrical engineering: 
the transverse oscillating electric field in the pulse is transformed into a longitudinal qua- 
sistatic electric field. The dimensionless amplitude of the longitudinal field is ay 150 
corresponding to E\\ = 4.8 x 10 14 V/m x (1 /im/A). The typical distance over which 
charge separation occurs is comparable with the initial thickness of the foil and is much 
smaller than the transverse size of the region that is being pushed. The proton layer is ac- 
celerated by this longitudinal field. We note that the laser pulse frequency in the reference 
frame comoving with the accelerated plasma region decreases as time progresses so that 
the accelerating foil become less transparent with time. 

As seen in the cross-section of the electric field component E z in Fig. 2, the thick- 
ness of coloured stripes, which corresponds to half of the radiation wave length, increases 
from left to right in the reflected part of the pulse (along the x-axis). This increase is 
weaker at the periphery (in the transverse direction). This 'nonuniform red shift' results 
from the Doppler effect when the laser light is reflected from the co-propagating rela- 
tivistic mirror which accelerates and deforms in time. The red shift testifies that the laser 
pulse does indeed lose its energy by accelerating the plasma mirror. In this stage, the 
foil is transformed into a "cocoon" inside which the laser pulse is almost confined. The 
accelerated protons form a nearly flat "plate" at the front of the "cocoon" as is seen in 
Fig. 3. 

Fig. 4 shows the maximum proton energy versus time. This dependence is initially 
linear and, at later times, the maximum proton energy scales as t 1//3 as predicted by the ID 
analytical model of Sec. 2.1. The protons in the plate structure are accelerated according 
to the RPD regime. These results provide numerical evidence of the fact that the RPD 
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Figure 2: Spatial distribution of the proton density and the z component of the electric 
field during the acceleration phase (see text for details). 




Figure 3: Spatial distribution of the proton density at t = 100. 
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Figure 4: Maximum proton energy versus time. The continuous line represents the ana- 
lytical expectation. 

acceleration mechanism should appear already for laser intensities of 10 23 W/cm 2 . The 
time evolution is hydrodinamically stable and the acceleration highly efficient. In com- 
parison with the experimental results of present day Petawatt lasers, this example predicts 
yet another astonishing advantage of the Exawatt lasers [16,30], besides those described 
inRefs. [16,31]. 

In Fig. 5 we show the proton energy distribution obtained in the 3D PIC simulations 
and their transverse emittance. The energy spectra have been calculated for the particles 
in the region near the beam axis ("plate") within a l/im radius. We see that the protons 
have a finite-width spectrum, localized within the interval 1.3 GeV <£ p kin< 3.2 GeV . 
The number of protons, integrated over energy, in the dashed region in Fig. 5 is equal 
to 2.7 x 10 10 particles per /im 2 . The proton transverse emittance is almost constant and 
equals ~ 1CT 2 7T mm mrad. 

The numerical studies depicted above are extremely challenging even for large par- 
allel computer facilities. In order to reduce complexity, the study has been carried out 
with laser pulse of relatively small focal spot. In addition, the dynamical evolution has 
been followed up until the t 1 ^ asymptotic behaviour is reached (i.e. before the complete 
laser-plasma decoupling). The overall laser energy to proton kinetic energy conversion 
efficiency (e) at that time is 40%. Extrapolation up to the time of decoupling indicates 
that an energy conversion efficiency of 57% can be reached. On the other hand, due to 
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Figure 5: Proton energy distribution (red) and transverse emittance (blue) near the beam 
axis at t = 80. 

the small focal spot, the RPD regime is operational only the central area: here a nearly- 
monochromatic spectrum is observed. The peripheral area ("cocoon") result in a nearly- 
thermal spectrum of protons. Hence, if we consider only the RPD accelerated protons 
at time of decoupling, the efficiency drops to 11%. This value is the most conservative 
estimate of conversion efficiency and it is not expected to hold for large focal spots S: the 
"cocoon" is the outcome of a border effect due to the finite waist of the pulse and will not 
scale linearly with S. Finally, note that the intensity and the duration of the pulse has not 
been optimized for the production of protons in the few GeV range but has been chosen to 
demonstrate the possibility of highly relativistic ion generation [17]. Awaiting for further 
numerical studies, in the following sections the proton and neutrino fluxes are provided 
as a function of the repetition rate and realistic expectations are discussed below. 

2.3 Proton fluxes 

A proof of principle of the RPD acceleration mechanism is at the borderline of current 
technology, but the possibility of using this technique to overcome the limitation of tra- 
ditional proton accelerators faces many additional difficulties. Present day systems based 
on Chirped Pulse Amplification [32] (CPA) are able to deliver intensities of the order of 
10 22 W/cm 2 . It can be expected that the intensity needed to trigger the RPD mechanism 
will be reached in the near future since for intensities up to 10 23 W/cm 2 the saturation flu- 
ence remains below the damage threshold and we can still profit of CPA for short-pulse 
generation. On the other hand, all present high power lasers operate at very low repetition 
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rate. This is a classical problem e.g. in inertial fusion: here, the basic principles could be 
demonstrated through the construction of optical cavities delivering up to 1.8 MJ as the 
National Ignition Facility (NIF) in US or the Laser Megajoule (LMJ) in France. How- 
ever, the ultimate use of ICF to produce electric power will require repetition rates of the 
order of tens of Hz, an increase of several order of magnitude compared to the shot rate 
achievable with state-of-the-art fusion laser technology. This rate cannot achieved with 
flashlamp-pumped neodymium-doped glass lasers that requires a significant interpulse 
cooling time. There is, presently, a very large effort to find alternative solutions [33]. We 
note here that the solution of the problem of thermal stability of the system for energy 
yield of tens of MJ/s would simultaneously provide an appropriate driver for ICF operat- 
ing with ~ 2MJ laser pulse at a repetition rate of 0( 10) Hz and, through the exploitation 
of CPA and the RPD acceleration mechanism described above, an unsurpassed proton 
acceleration facility for hadron and spallation neutron production operating at energies of 
~1 GeV (laser pulses of 16 kJ/e) and repetition rates of the order of a few kHz. As a 
purpose of illustration, Tab. 1 compares the main parameters of present and future pro- 
ton drivers at the GeV energy range. The LAMPF and ISIS beam dump facilities have 
given neutrino sources for the LSND [34] and KARMEN [35] experiments. Currently, 
two projects (JAERI in Japan [36] and SNS in USA [37]) are under construction, while 
the European Spallation Source project (ESS) [38] and the CERN Superconducting Pro- 
ton Linac (SPL) [39] are still pending approval. A RPD laser-driven acceleration facility 
(LAF) operating at 1 kHz and providing 10 14 proton/pulse (i.e. with a "plate" radius of 
~34 fim) is also shown. It corresponds to a energy yield of about 50 MJ/s at e = 0.3. It is 
a remarkable fact that the RPD acceleration mechanism allow a close synergy between the 
technological needs for a ICF laser driver and an ultra-intense proton driver for nuclear 
and particle physics applications. 

3 Neutrino beams from n + and fi + decay-at-rest (DAR) and decay-in-flight (DIF) 
3.1 The 7r + and fi + decay chains 

In a neutrino beam-line based on a dump of low energy protons into a passive material, 
neutrinos arise from both pion and muon decays. The production of kaons or heavier 
mesons is negligible if the proton energy is sufficiently low (E p < 3 GeV). Therefore, the 
neutrino beam does not suffer from v e and v e contaminations due to kaon decays. The 
pion decay modes are 7r + — > /x + z/ M , 7r + — > e + u e , ir~ — > /i~z> M and ir~ — > e~v e but the 
decays into electrons are strongly suppressed. The muon decay modes are fi + — > e + z/ e z/ M 
and fi~ — > e~z/ e z/ M . Almost all /x+ stop before decaying and produce a Michel spectrum 
for u e and while the fx~ are captured in orbit. The ir + decay occurs both with the pion 
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Table 1 : Main parameters of present and future proton drivers at the GeV energy range. 
Pot = proton on target. 





Energy 


pot/pulse 


Beam 


Repetition 


Pulse 


Target 




(GeV) 


xlO 13 


current 


rate 


width 




ISIS 


0.8 


2.5 


200 //A 


50 Hz 


100 ns 


Tantalum 


LAMPF 


0.8 


5.2 


1 mA 


120 Hz 


600 /is 


Water/High Z 


JAERI 


3.0 


8.3 


333 /lA 


25 Hz 


1 /IS 


Hg 


SPL 


2.2 


15 


1.8 mA 


75 Hz 


2.2 ms 


various 


SNS 


1.0 


15 


1.4 mA 


60 Hz 


695 ns 


Hg 


ESS 


1.3 


84 


6.7 mA 


50 Hz 


1.4 /is 


Hg 


LAF 


1.0 


10 


16 mA 


1 kHz 


<1 ps 


Water/High Z 



at rest, providing a mono-energetic neutrino spectrum, and in flight. The ratio DAR/DIF 
depends on the material and the geometry of the target. Proton-rich targets (e.g. water) are 
employed to obtain large pion yields. Early stopping of the mesons is achieved positioning 
a dense dump just after the water vessel. The distance between the water target and the 
stopper can be tuned to optimize the DAR/DIF ratio. 

3.2 Neutrino energy spectra and beam composition 

In order to estimate the expected neutrino energy spectra and beam composition, we re- 
fer to the setup of the LSND experiment, operated at LAMPF from 1993 to 1998 with 
800 MeV protons impinging into a water target followed by a copper beam stopper. The 
precise evaluation for a laser driven facility should include the secondary yield of the 
nearly monochromatic distribution plus the higher energy tail (Fig. 5) and the contribu- 
tion of the quasi-thermal spectrum. The former will be peaked in the 1-2 GeV range to 
make fully operative the RPD mechanism. The size of the latter will depend on the co- 
coon/plate ratio, as described in Sec. 2.2. Clearly, a precise determination of the fluxes 
is beyond the scope of this paper; however, we note that the use of LAMPF data implies 
a significant underestimation (~ 50%) of the n + yield and a small underestimation of 
the background from DIF which can be reduced through a dedicated optimization of the 
target-stopper distance. In the following we refer to a laser-driven facility providing 10 14 
protons-on-target (pot) per pulse at different repetition rates with a nearly monochromatic 
spectrum corresponding to the LSND setup 4 . 

4 Full flux calculations are available in [40]. In fact, all proposed future neutron spallation sources have 
as target material liquid mercury. The reason for such a choice is twofold: mercury is liquid at room 
temperature. Therefore, its recirculation allows a more efficient power dissipation with respect to solid 
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In the LSND target configuration the DAR/DIF ratio turned out to be 97%/3% [41- 
43]. The leading decay chain ir + — > — > e + z/ e z/ M z/ M does not contain 9 e and offers 
a unique opportunity to test the occurrence of — > z/ e transitions. The decay of n~ 
might lead, in principle, to a large z/ e contamination. However, three factors contribute 
to its suppression. At these proton energies, the ir + production rate is larger than ir~ 
by about a factor 8. Moreover, negative pions which come to rest are captured before 
they decay: in fact, at LAMPF, only 5% decay in orbit and, hence, contribute to the u e 
background. Finally, almost all negative muons arising from the decays in flight come to 
rest in the beam dump before decaying; most then undergo the reaction fi~N u^e~ 
that leads to with energy below 90 MeV, leaving only 12% of them to decay into v e . 
Hence, the relative v e yield of the ix~ decays at LSND, compared to the n + decays, is 
~ (1/8) x 0.05 x 0.12 pa 7.5 x 10~ 4 . Clearly, the level of v e contamination is much 
smaller than the intrinsic v e contamination (a few %) of a high energy beam from % 
decay (e.g. the Superbeams). 

The 7r + decay chain provides an intense source of monochromatic muon neutrinos 
(E = 29.8 MeV). Due to the short n + lifetime (r n ~ 26 ns) these neutrinos closely follow 
the beam time profile. This fact opens up the possibility to detect — > v e oscillations: 
the — > v e oscillations can be temporally separated from the events due to /i + DAR 
electron neutrinos, which appear on a time scale of few /is due to the muon lifetime 
(t m ~ 2.2 /is). In order to exploit the well defined time structure of a DAR beam, it is 
mandatory to have a proton pulse width comparable with the pion life-time. For a laser- 
driven facility this requirement is easily fulfilled, the proton pulse temporal spread being 
of the order of 1 ps. 

Finally, the shape of the neutrino flux from ir + and /x + decay at rest (DAR) is well 
known and it is shown in Fig. 6. Therefore, only the absolute amplitude has to be deter- 
mined from experiments and simulation. 

The simulation of the expected neutrino fluxes, both at LAMPF and ISIS, was per- 
formed by using the pion yield in proton-target interaction measured in a dedicated ex- 
periment [44]. For details on the neutrino flux simulation and the associated uncertainties 
we refer to [45,46]. The main source of systematic error is associated with the pion 
yield in proton-target interaction and it has been estimated to be about 6%. This has to 
be compared with the total systematic error associated to the neutrino flux from decay at 
rest that has been estimated to be of about 7% [45]. For neutrinos from decay in flight the 

target. Moreover, the high atomic number gives a source of numerous neutrons. For neutrino applications 
only the first motivation holds. Indeed, high-Z elements are less efficient in producing neutrinos that low-Z 
elements. On the other hands, high-Z elements allow a strong suppression of DIF neutrinos. In water-based 
targets this suppression, if needed, can be partially recovered modifying the distance between the water 
vessel and the stopper. 
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Figure 6: Neutrino energy spectra from n + — > yU + z/ M — > e + z/ M z/ e z/^ decay at rest, 
systematic error has been estimated to be about 15% [42]. 
4 The neutrino oscillation channels and the detector 

In principle neutrinos from DAR and DIF decays could allow simultaneous investigation 
of 

• — > v e oscillations through the identification of DAR transitions ("DAR anal- 
ysis"); 

• — > v e oscillations from high energy DIF ("DIF analysis"); 

• v p — > v e oscillations from DAR temporally separated with respect to DAR v e 
("time analysis"). 

So far these channels have been used to explore neutrino oscillations at a Am 2 of about 
1 eV 2 (E ~ 10 1 MeV and L ~ 10 1 m). The LSND experiment gave a positive result 
(there is a claim for an excess of events with an electron in the final state induced by 
(^/i) oscillations into u e (z/ e ) [41-43]), while the KARMEN experiment gave a nega- 
tive result [47]. The need for a check of these experiments caused the proposal of new 
projects like MiniBooNE [48] (currently data taking) and of new experiments at the Neu- 
tron Spallation Source [49]. However, there are no proposals to search with DAR and/or 
DIF neutrinos for oscillations at a baseline of few kilometers to test sub-dominant v e and 
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V e appearance modes at the atmospheric scale. This test is unfeasible with present accel- 
erators. Much higher intensities ((9(20 mA) or more) are needed to overcome the large 
suppression due to the smallness of the v e CC cross-section at these energies. As noted 
before, a laser driven facility could offer this opportunity. 

The optimization of the detector coupled with this facility is beyond the scope of 
this paper. Here, we will consider a detector with a technology similar to LSND and with 
a fiducial mass comparable with Super-Kamiokande. This corresponds to 1.3 x 10 33 free 
protons (17 kton of CH 2 ). Similar detectors based on liquid scintillator have been recently 
proposed for low-energy neutrino astronomy and proton decay [50]. In the following, 
when appropriate, we refer to the LSND experimental and Monte Carlo studies to assess 
the physics performance of the apparatus [41-43]. 

The LSND detector consisted of a cylindrical tank filled with liquid scintillator. The 
composition was chosen to be sensitive to both Cerenkov light from electrons and rela- 
tivistic muons and scintillation light from all charged particles. The light was detected 
through photomultiplier tubes (PMT's) covering 25% of the detector surface. PMT time 
and pulse-height signals were used to reconstruct the track. The Cerenkov cone for rela- 
tivistic particles and the time distribution of the light, which is broader for non relativis- 
tic particles, gave excellent separation between electrons and particles below Cerenkov 
threshold. 

5 The Pfj, — > v e decay at rest analysis 

The search for — > v e oscillations is performed by using from fi + decay at rest. 9 e are 
detected through the reaction v e p — > e + n (plus a small contamination of u e C — > e + B n) 
followed by the neutron capture reaction np — > dj, the 7 energy being 2.2 MeV. A can- 
didate P e events consists in one identified electron with energy in the 20 < E e < 60 MeV 
range 5 and one associated gamma. The electron identification efficiency at LSND is 42% 
with a relative systematic error of 7%. Note that the e + inefficiency is dominated by the 
need of vetoing Michel electrons from cosmic muons; in particular the dead-time of the 
veto accounts for an electron efficiency reduction of 24% [43]. The size of this back- 
ground strongly depends on the shallow depth of the experimental area where LSND is 
located and would be significantly suppressed at deeper locations. The correlated photon 
is identified combining the information from the number of PMT's hits associated with 
the 7, its distance from the positron and the time interval between the e + and the 7 which 
exploits the 186 /is delay for the neutron capture in mineral oil. The LSND analysis is 
based on a likelihood function i? 7 whose i? 7 > 10 cut corresponds to an efficiency of 

5 /3 decays of cosmogenic 12 B prevent the use of the candidate with E < 20 MeV. 
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Figure 7: Energy spectrum of the fully oscillated — > v e (dotted line) before the selec- 
tion cuts. The two-family oscillation probability for Am 2 = 2.5 x 10~ 3 eV 2 at L = 11 km 
(continuous line) and the corresponding convoluted — > v e spectrum (dashed line) are 
also shown. 

39% and a contamination due to accidentals at the level of 0.26%. The corresponding 
relative systematic uncertainty does not exceed 7%. Fig. 7 shows the energy spectrum of 
— ► ^ e oscillated neutrinos for 100% conversion probability before any selection cut: 

$^(E) a DeP ^ e+n (E) (23) 

i.e. the product of the v^E) flux from /i + DAR and the u e p — > e + n CC cross-section. 
The two-family oscillation probability for Am 2 = 2.5 x 10~ 3 eV 2 at L = 11 km and the 
corresponding 

a DeP ^ e+n {E)P(P^ -> P e )(£) (24) 

product is also shown. The corresponding cross- section, weighted with the DAR spec- 
trum, is 0.95 x 1(T 40 cm 2 . 
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Contaminations to the v e sample arise from beam-related interactions with neutrons 
in the final state, beam events without neutrons and background produced by cosmics and 
radioactivity ("beam unrelated"). 

The main beam-related backgrounds with final state neutrons is v e p — > e + n from 
fx~ decays at rest. As noted before, the fx~ DAR yield is suppressed by a factor 7.5 x 1CT 4 
w.r.t. /x + DAR. The v e p — > e + n cross-section weighted according to the P e DAR spec- 
trum is 0.72 x 10~ 40 cm 2 [51]. The fraction of events with energy greater than 20 MeV is 
0.806 while the electron and 7 efficiencies are the same as for the signal sample. The sec- 
ond most important source of beam-related background events with correlated neutrons 
is the misidentification as P e events of p — > /i + n CC interactions from 7r _ DIF. Be- 
cause of the energy needed to produce a such a must arise from a decay in flight. 
The final state considered are v^p — > ji + n or (less often) v^C — > ji + nX, followed by 
/i + — > e + v e i / fl . In most of the cases, the muon is missed because it decays at very large 
times compared with 7^ ~ 2.2 /is or the deposited energy is below the phototube thresh- 
old. The latter can occur either because the muon is too low in energy or it is produced 
behind the phototube surfaces. For this background, the flux weighted cross-section is 
4.9 x 10~ 40 cm 2 , the fraction of muons in the tail of the lifetime distribution (r > 12 /xs 
or with very low kinetic energy (T < 3 MeV)) is 2.6%. The positron efficiency is 42%, 
the fraction of events with E > 20 MeV is 81.6% and, again, the efficiency for correlated 
7 is 39%. Other source of misidentification are muon decays in the tail of the lifetime dis- 
tribution, prompt decays to electrons so that the /i and the e are collected in a single event 
and muon lost by trigger inefficiencies. The whole background from /j, misidentification 
has been computed according to the results of [43]. However, the DIF flux has been cor- 
rected keeping into account the larger distance of the detector (L = 11 km versus 30 m): 
at large baselines the DIF fluxes are similar to the corresponding fluxes at the center of 
LSND corrected for the L~ 2 suppression term. 

The main source of beam related background without correlated neutrons is v e 12 C — 
e~X scattering. The corresponding average cross-section is 1.5 x 10~ 41 cm 2 [52]. For 
an electron reconstruction efficiency of 0.36 [41], the fraction of events with E > 20 
MeV is 46% while the electron efficiency and the accidental 7 efficiency at LSND are, 
respectively, 42% and 0.26%. Other sub-dominant sources are discussed in [41]. 

Beam unrelated background results mainly from unvetoed cosmic interactions in 
delayed coincidence with accidental photons. It strongly depends on the veto quality, the 
depth of the detector and the rate of accidental photons. Moreover, it can be suppressed 
if the duty cycle of the beam is sufficiently low. For the case of LSND, the poor time 
structure of the beam and the shallow depth of the detector does not allow an effective 
suppression of this background, which is estimated during the beam-off data taking and, 
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hence, subtracted. On the other hand, the time structure of ISIS allows a suppression 
of more than two order of magnitude compared with LSND in the — > 9 e channel, so 
that KARMEN does not suffer from this contamination. For the laser-driven facilities 
considered here the time structure is very well defined (see Sec. 2 and [18]) and the 
beam-off background could be non negligible (at the depth of LSND) only for very high 
repetition rates (~ 10 kHz). At larger depths, this contamination is negligible for any 
realistic repetition rate and, hence, it is not considered in the present analysis. 

6 Sensitivity to 6>i 3 

In order to get information about the magnitude and phase of the U e3 term of the PMNS 
mixing matrix, terrestrial experiments explore sub-dominant effects in the neutrino tran- 
sition probabilities at the atmospheric scale which, in general, are suppressed by at least 
one power of a = Am^/IAm^l 6 . Since matter effects are negligible for baselines of 
the order of L = 10 km, the — > v e oscillation probability can be expressed as [53]: 

P{ v n v e) — sin 2 26*13 sin 2 6*23 sin 2 A 
—a sin 2#i3 sin 5 cos #13 sin 2#i2 sin 26*23 sin 3 A 
—a sin 26 ) 13 cos 5 cos #13 sin 26 , 12 sin 26* 2 3 cos A sin 2 A 

+a 2 cos 2 6* 23 sin 2 26 12 sin 2 A 

= Oi + 2 + 3 + 4 , (25) 

A being the oscillation phase Am| 3 L/AE U in natural units (c — k— 1). If the energy and 
the baseline of the experiment is chosen to fulfill A ~ n/2 (oscillation maximum), the 3 
term is suppressed. Similarly, the term 4 can be neglected unless sin 2 26 13 ~ (9(10~ 4 ). 
Hence, the — > v e transition probability can be expressed in the simplified form: 

Piv^ -> u e ) ~ sin 2 26»i3 sin 2 6» 2 3 sin 2 A 

— a sin 26*13 sin 5 cos 9i 3 sin 26*12 sin 26*23 sin 3 A (26) 

Its CP conjugate — > v e is therefore 

P(y ^ -»• v e ) ~ sin 2 26 l3 sin 2 6*23 sin 2 A 

+ a sin 26*13 sin 5 cos 6*i 3 sin 20 12 sin 26* 2 3 sin 3 A (27) 

and a simultaneous measurement of u e and u e gives access 7 to the magnitude and phase 
of the Ues entry. 

6 A global analysis of all available neutrino oscillation data gives for a a best fit value of 0.026, while 
the 3a allowed range is 0.018 < a < 0.053 [6]. 

7 This determination is not unique due to the S — > n — S and 6*23 — > n/2 — 623 invariance of the transition 
probabilities [54,55]. 
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The analysis described in Section 5 allows a determination of Pip^ — > P e )- To 
ease comparison with other proposed facilities, we interpret Piy^ — > z> e ) assuming no CP 
violation in the leptonic sector (5 = 0), # 2 3 = vr/4 and Am| 3 = 2.5 x 1CT 3 eV 2 , so that 



P(i/ (U — > u e ) ~ sin 2 2^13 sin 2 #23 sin 2 A 



^ sin 2 2^13 sin 2 



1.27 



Am 2 3 (eV 2 ) L(km) 
E(GeV) 



(28) 



and performing a two parameter fit of Am 2 3 and sin 2 20 13 . The fluxes of from n + DAR 
are computed for a facility providing 10 14 pot/pulse with variable repetition rates. The ir + 
yield considered here is 0.09 7r + /pot, corresponding to the LSND setup [45]. As discussed 
in Sec. 3.2, this estimate is rather conservative for a laser-based facility. In 5 years of data 
taking, assuming 6 months of operation and 50% beam on time, the overall integrated flux 
at 11 km is 2.3 x 10 7 x R u/cm 2 , R being the repetition rate in Hz. At R = 1 kHz, we 
expect 473 P e CC events assuming 100% — > v e conversion rate. The corresponding 
background is below 0.6 events (0.28 from pT DAR, 0.16 from 7r DIF with the muon 
misidentified and 0.16 from events without correlated 7). The systematic errors for signal 
and background are inferred from the LSND analysis and shown in Table 2. Fig. 8 shows 
the 90% confidence level (CL) exclusion limit for sin 2 2#i3 in the occurrence of the null 
hypothesis (9 13 = 0) as a function of the repetition rate. The corresponding sensitivities 
coming from the — > v e appearance search for MINOS [56], CNGS [57], JPARC- 
SK [9] and JPARC to Hyper-Kamiokande [9] are also shown, together with the present 
CHOOZ [2] limit from v e disappearance. The two-parameter exclusion region at 90% CL 
for a 10 kHz facility is depicted in Fig. 9. 



7 h>^—>i> e appearance searches 

In LSND — * v e oscillations were searched for by exploiting the flux from decay in 
flight. As discussed in Section 3, about 3% of the n + decay in flight originating a beam 
with average energy of about 90 MeV. v e are detected through the reaction v e e — > e X 
and requiring the electron to have energy in the range 60 -j- 200 MeV. The corresponding 
value of the cross-section for < E v >~ 9 MeV is about 15 x 10~ 40 cm 2 [52], i.e. a 
factor 15 larger than v e inverse /3-decay (see Section 5). Taking into account the different 
detection efficiencies [42] we expect a ratio of fully oscillated DIF/DAR events of about 
one third with a signal to noise ratio comparable to both analyzes. Note, however, that the 
energy of DIF neutrinos does not match the maximum of the oscillation probability for 
the baseline considered (L — 11 km for Am| 3 = 2.5 x 10 -3 eV 2 ). 
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Figure 8: Exclusion limit at 90% CL for sin 2 26>i 3 (two parameter fit) as a function of 
the repetition rate. The corresponding integrated protons on target for 5y of data taking 
are also shown. The horizontal bands indicate the corresponding limit from — > v e ap- 
pearance search at MINOS, CNGS, JPARC to Super- Kamiokande and JPARC to Hyper- 
Kamiokande; the upper band shows the present CHOOZ limit from v e disappearance 
search. 
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Table 2: Breakdown of systematic errors for the — > v e DAR analysis. 



Channel 


Source 


Error (%) 


signal 


flux 


7 




cross-sec 


_ 




eff 


7 




7 id 


7 


v e from DAR 


flux 


15 




cross-sec 






eff 


7 




7 id 


7 


from DIF ir~ 


flux 


15 


eff + cross-sec 


41 




7 id 


7 


u e 12 C from DAR ji + 


inclusive 


17 




Figure 9: Exclusion limit at 90% CL in the (Ara| 3 , sin 2 26*13) plane for 10 kHz repetition 
rate. 
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The time analysis searches for — > v e oscillations by looking for mono-energetic 
v e (E u = 29.8 MeV) in a short time window after the proton pulse. The length of the time 
window is determined by the pion life-time (r^ = 26 ns). The v e are detected through the 
CC reaction of u e onto 12 C giving rise to an electron and a 12 N (ground state) followed 
by the /3-decay 12 A^. S . — > 12 C e + v e with 15.9 ms decay time. This analysis was found 
particularly appealing for the following reasons: the main background is induced by v e 
from fast yU + -decay within the time window, but can be precisely measured and subtracted; 
due to v e 12 C — > e~ 12 N g , s , interactions from fi + decays outside the time window, the full 
oscillation expectation is normalized by the experiment itself. However, the sensitivity 
of this channel is essentially limited by the small cross-section of the involved reaction 
(4.95 x 10~ 42 cm 2 , at E v = 29.8 MeV, to be compared with the inverse /3-decay reaction 
whose cross-section is 0.95 x 10~ 40 cm 2 , convoluted over the whole DAR spectrum). 
Summarizing the search for z/ M — > v e oscillations is very interesting for CP-violation 
studies, but it is not at the peak of the oscillation probability in the DIF analysis (E/L ~ 
0.1 GeV/1 1 km ~ 10~ 2 ; A ~ 0.34) and it is limited by statistics in the time analysis. The 
impact of this channel on CP-violation studies is currently under investigation and will be 
the subject of a forthcoming paper. 

8 Conclusions 

The current debate for multipurpose facilities aimed at high precision studies of neutrino 
oscillations drove the authors towards the study of non-conventional neutrino sources. In 
particular, we noted that the possibility to accelerate efficiently protons in the GeV energy 
range through relativistic laser-plasma interactions opens up interesting opportunities for 
the development of a new generation of proton drivers. In this paper we discussed a 
radiation pressure dominated (RPD) mechanism for relativistic proton acceleration. This 
mechanism is highly efficient compared to previous proposals and could allow a close 
synergy between present R&D finalised to energy production through inertial confined 
fusion and the wealth of applications related to high intensity multi-GeV drivers (see 
Appendix). Moreover, we demonstrated that this facility could allow for the first time 
the study of subdominant — > v e oscillations at the atmospheric scale with neutrinos 
produced by n + decays at rest or in flight. 
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Appendix 

In the last decade, interest towards the construction of high intensity proton drivers in the 
few GeV range has steadily grown. The main reason is connected with the wide range 
of applications that can be simultaneously accessed by these facilities. The most intense 
neutron beams (spallation sources) are currently produced by bombarding mercury targets 
with energetic protons from a large few-GeV proton accelerator complex. As it is well 
known, applications range from chemistry to crystalline and disordered materials studies, 
superconductivity, polymers and structural biology investigations. A discussion of the 
physics case for intense spallation sources can be found in [58]. A wide physics program 
is accessible in neutrino physics, beyond the oscillation issues discussed above [49]. The 
availability of a high power, high duty factor proton beam could provide opportunities in 
stopped muon physics, in search for rare decays as /i — > cy, /i — > eee or muon conversion 
/iN — > eN, and improvements in muon decay properties. The neutrino energy range is 
appropriate for neutrino-nucleus cross section measurements of relevance to supernova 
astrophysics: dynamics, nucleosynthesis and terrestrial supernova v detection. Similarly, 
measurements of neutrino-nucleus cross sections open the possibility to study interesting 
nuclear structure issues related to the weak interaction as the ratio of the axial to vector 
coupling constants and the search for non-standard contributions. Finally it's worth noting 
that oscillation studies at the Am 2 ~ 1 eV 2 scale will become a major priority in v physics 
in case of confirmation of the LSND v e appearance claim. 
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